LONG-TIME ASYMPTOTICS FOR FULLY NONLINEAR 
HOMOGENEOUS PARABOLIC EQUATIONS 

SCOTT N. ARMSTRONG AND MAXIM TROKHIMTCHOUK 

Abstract. We study the long-time asymptotics of solutions of the uniformly 
parabolic equation 

u t + F(D 2 u) = in 1" x 1 + , 

for a positively homogeneous operator F, subject to the initial condition 
u(x, 0) = g(x), under the assumption that g does not change sign and possesses 
sufficient decay at infinity. We prove the existence of a unique positive solution 
and negative solution <J>~, which satisfy the self-similarity relations 

$±0,*) = A a± * ± (A 1 / 2 a;, Xt). 

We prove that the rescaled limit of the solution of the Cauchy problem with 
nonncgative (nonpositive) initial data converges to <J> + (<t> — ) locally uniformly 
in R n X R_|_. The anomalous exponents »+ and ct~ are identified as the 
principal half-eigenvalues of a certain elliptic operator associated to F in R n . 



1. Introduction and main results 

The connection between the scaling invariance of the mathematical expressions 
for certain physical laws and the asymptotic behavior of physical phenomena is of 
fundamental importance to the study of mechanics. In this work, we have in mind 
the study of self-similar solutions of diffusion equations and their relationship to 
the long-time behavior of solutions to the Cauchy problem. A classical example is 
the heat equation 

(1.1) u t ~Au = in I" x I + , 

which is invariant with respect to any scaling (x, t) t— > (a 1 ^ 2 x^ at) , with a > 0. It is 
well-known that a solution u(x, t) of equation (|1.1[) with nonnegative and integrable 
initial data will converge in a certain sense to a multiple of the Gaussian kernel 

$(x,t) := (47rf) 2 e « . 

Precisely, the rescaled solutions u a (x,t) := a n l 2 u {o 1 l 2 x,o£) will converge locally 
uniformly to C$>{x, t), where the constant C is given by 

C = I u(x, 0) dx. 
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The Gaussian kernel <& satisfies the relation 

and for this reason $ is called a self-similar solution of 

Another important example of a diffusion equation is the generalized porous 
medium equation 

(1.2) u t -A*(u)=0 in M" x R +! 

where 'J is an increasing function of u. The study of self-similar solutions of (jl.2p 
and the long-time asymptotics of solutions to the Cauchy problem is well-developed, 
particularly in the case that 'J/ (it) = u m for m > 1. It originated with the celebrated 
work of Barenblatt, Zel'dovich, and Kompaneets, and was later developed Fried- 
man, Kamin, Vazquez, and others (see Vazquez [24| for a well-written introduction 
to the subject as well as a comprehensive list of references). 

While the literature on self-similar solutions and asymptotics of diffusion equa- 
tions is vast, relatively little is known in the case of the fully nonlinear parabolic 
equation 

(1.3) u t + F(D 2 u) = in R n x R + . 

Here F is a positively homogeneous, uniformly elliptic operator. Important exam- 
ples of (jl.3p include the parabolic Bellman and parabolic Isaacs equations, which 
arise in the theory of stochastic optimal control and stochastic differential games, 
respectively. In the article, we show that equation (|1.3[) possesses a unique positive 
and negative self-similar solution, and that these self-similar solutions characterize 
the long-time asymptotic behavior of solutions of the Cauchy problem with initial 
data which does not change sign. 

Theorem 1.1. Assume that F satisfies (I2.2[) and (|2.3[) . below. Then there exist 
unique constants a + = a + (F) > and aT — a~(F) > 0, for which the uniformly 
parabolic equation (|1.3|) possesses solutions $ + > and $~ < 0, satisfying the 
relations 

<S> + (x,t) = ct q+ $+ (V /2 .t, eri) and <f>~(x, t) = a a ~ (a 1/2 x, at) 

for all a > 0, and such that for some constants C, a > 0, 

$ + (a;,l), -$ _ (a;, 1) < C*cxp(-a|x| 2 ) for all x G M™. 

Moreover, the solutions $ + and $~ are unique up to multiplication by a positive 
constant. 

We call the numbers a + and a~ the positive and negative anomalous exponents 
of the operator F, respectively. In contrast to the case of a linear parabolic equation, 
a 7^ n/2 and a + ^ or , in general. We will see in Section [3] that they are the 
principal half-eigenvalues of an elliptic operator in R n . The functions < I >+ (-, 1) and 
$~(-, 1) are the corresponding principal half-eigenfunctions. 

Our second main result characterizes the long-time behavior of solutions of equa- 
tion (|1.3p with initial data which does not change sign and exhibits sufficient decay 
at infinity. 
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Theorem 1.2. Let the hypotheses of Theorem li.il be in force, and consider a 
viscosity solution u E C(R™ x R+) of the uniformly parabolic equation (|1.3|) such 
that |u(ir,0)| < Co exp(— B|cc| 2 ) for some constants B,Cq > 0. Ifu(-,0) > and 
u(-,0) ^ 0, then there exists a constant C* > such that the rescaled solutions 
given by 

u a (x, t) := a a+ u (a 1/2 x, at) , (x, t) E R" x R + 

converge locally uniformly in R™ x R + as a — > oo to the function C*& + . Likewise, if 
u(-, 0) < 0, u(-, 0) ^ 0, then there exists a constant C* > such that the functions 

u a (x,t) := a a ~u (o- 1/2 a;, cri) , (x,t) E R" x R + 

converge locally uniformly to C*<fr~ as a — > oo. 

The Barenblatt equation of elasto-plastic filtration 
(1.4) Mj-max , ^=0 in R" x R+, 



1 - 7 1 + 7 . 

where < 7 < 1, is a particular example of a fully nonlinear parabolic equation of 
the form (|1.3[) . It arises in filtration theory as a model for an elastic fluid flowing 
through an irreversibly deformable elasto-plastic porous medium (see Barenblatt, 
Entov, and Ryzhik Q). Kamin, Peletier and Vazquez [T3| proved Theorems 11.11 and 
1.21 for equation (|1.4|) . using different methods from those employed in this paper. 
Their argument makes use of the rotational invariance of the equation to reduce the 
problem to an ODE. The authors then employ a shooting argument in the phase 
plane to demonstrate the existence of self-similar solutions of equation (| 1 . 4|) . Their 
proof of the asymptotic convergence of rescaled solutions of the Cauchy problem to 
a multiple of the self-similar solution relies on a careful analysis of the self-similar 
solution at infinity. These estimates also rely on ODE theory in a crucial way, and 
in particular some clever applications of l'Hopital's rule. 



Our proof of Thcorcm l 1 . 2l cmplovs the "four-step method" of Kamin and Vazquez 
flit in a new setting. The idea is to squeeze the rescaled functions u° between two 
multiples of the self-similar solution and to show that the gap between them 



must vanish asymptotically. Similar to 171 ]. the primary difficulty in carrying out 



this analysis is controlling the behavior of the solutions as t 1 / 2 |x| — > 00. This is 
overcome by the construction of a special comparison function. 

Self-similar solutions of equation (|1.4[) have also been considered by Caffarelli and 
Stefanelli @], and the dependence of the anomalous exponents on the parameter 7 
has been examined by Goldenfeld, Martin, Oono, and Liu 12] as well as Aronson 
and Vazquez Q. We also wish to mention the interesting work of Hulshof and 
Vazquez [15| , who studied the long-time asymptotic behavior of viscosity solutions 
of an equation similar to (II. 4|) . but with A(u m ) replacing both instances of Au. 

As we were completing a final revision of this paper for publication, we became 
aware of a new preprint by Meneses and Quaas [201 ] in which the existence assertions 
of Theorem 1 1.1 1 are also obtained. 

In Section [2] we state our hypotheses and recall the definition of a viscosity 
solution. Section [3] contains the proof of Theorem 11.11 an d in Section |4] we prove 
Theorem H~2l 
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2. Notation and Hypotheses 

Throughout this paper, we denote (0,oo) by R+. The set of n-by-n real sym- 
metric matrices is S". For MeS" and < A < A, define the operators 

Vt A (M):= sup [- trace(AM)] and AM) := inf [- trace (AM)] , 
Ae[A,A] ' Ae [a,a] 

where [A, A] C S" is the set of positive definite matrices with eigenvalues contained 
in the interval [A, A]. The nonlinear operators "P^" A and » are called the Pucci 
maximal and minimal operators, respectively. For ease of notation, we will often 
drop the subscripts and write V + and V~ . A convenient way to write the Pucci 
extremal operators is 

(2.1) V+(M) = -A £ H - A £ and = -A £ ^ - A £ 

Aij>0 fj,j<0 iij>0 fj.j<0 

where /ii, . . . , fi n are the eigenvalues of M. 

We require our operator F : S" — > K to be uniformly elliptic in the sense that 
there exist constants < A < A, such that 

(2.2) V^ A (M — N) < F(M) - F(N) < V^ A (M - N) for all M, N G S n . 
We also require F to be positively homogeneous of order one: 

(2.3) F(r)M) = rjF(M) for all M G S", 7? > 0. 

Remark 2.1. If .F satisfies hypotheses (j2~2j) and |2~3)) . then the operator F(M) := 
— F(— M) satisfies ()2.2j) and (|2.3j) as well. This observation will simplify the proofs 
of our main results. For example, the existence of the negative self-similar solution 
for F can be deduced from the existence of the positive self-similar solution for 
F, and it is clear that a~(F) — a + (F). Likewise, to prove Theorem II .2\ we need 
only to assume nonnegative initial data and show convergence to a multiple of the 
positive self-similar solution, as the last statement follows from the former applied 
to F. 

Every differential equation and differential inequality in this paper is assumed 
to be satisfied in the viscosity sense. An introduction to the theory of viscosity 
solutions can be found in Crandall, Ishii, and Lions Q . For the convenience of the 
reader, we now recall the definition of a viscosity solution. 

Definition 2.2. Let G be an elliptic nonlinear operator, and £1 an open set in 
MP x R_|_. A function u G C(f2) is a viscosity subsolution (supersolution) of the 
parabolic equation 

(2.4) u t + G(D 2 u,Du,u,x) = in ft 
if whenever xq € ft and (p G C 2 (ft) are such that 

x i ► u(x) — ip(x) has a local maximum (minimum) at xq, 

we have 

u t + G(D 2 ip(x Q ),Dip(x ),u(xo) 1 x ) < (>) 0. 
We say that u is a viscosity solution of (|2.4p if it is both a viscosity subsolution 
and supersolution of (|2.4ll . Viscosity (sub/super)solutions of elliptic equations are 
defined analogously. The precise meaning of the differential inequality 

u t + G(D 2 u,Du,u,x) < (>) in ft 



ASYMPTOTICS FOR FULLY NONLINEAR PARABOLIC EQUATIONS 



is that u is a viscosity subsolution (supersolution) of (|2.4j) . 

Several times in this article, we will make use of the strong maximum principle 
and estimates for viscosity solutions of uniformly parabolic equations. For these 



facts we refer to Wang [25[ and Crandall, Kocan, and Swi§ch [10]. Analogous results 
for uniformly elliptic equations can be found, for example, in Caffarelli and Cabre 
0, Trudinger and Winter 

3. Existence of self-similar solutions 



In this section we prove Theorem ll.il To motivate our approach, suppose that 
there exists a solution $ > of equation (|1.3|) and an exponent a > which satisfy 
the relation 

(3.1) $0, t) = cr"$ (a 1/2 x, at) for every (x, t) € Q, a > 0. 
For t = 1 this reads 

(3.2) $0,1) = er"* (a 1/2 x,a) . 



Formally differentiate (|3.2p with respect to a at <r = 1 to discover that 

(3.3) $ t (a:, 1) + a$(x, 1) + -a; • £>$(£, 1) = 0. 

Inserting (|3.3|) into the PDE (|1.3p and rearranging, we derive the equation 
F (L» 2 $(x, 1)) - -x ■ D$(x, 1) = a$(x, 1) for all x e M". 

Defining </?(?/) := <&(?/, 1), we see that the pair (a, y>) is a solution of the elliptic 
eigenvalue problem 

(3.4) F (D 2 tp) — -y ■ Dip = aip in R". 

To prove Theorem ll.il we will reverse the calculation above. Studying the elliptic 
eigenvalue problem (|3.4p . we will show that there exists a > and <p > which 
satisfy (|3.4[) and such that ip decays at a suitable rate as \y\ — > oo. We will then 
define $ by 



$(x,<) :=r a cp[t 



-1/2 



and check that $ is a solution of (|1.3|) which also satisfies (|3.1I) . 

Another way to derive (13. 4p is to consider the continuous rescaling of (|1.3|) given 
by the change of variables 

y = t- 1/2 x, « = logt, ^(y,») = t a $(x,t), 

and then to look for a stationary solution (p(y,s) — tp{y). Thus the appearance 
of (|3.4p in our context is similar to the emergence of the so-called Fokker-Planck 
equation in the study of self-similar solutions of the porous medium equation (see 
Vazquez 0, Section 18.4]). An elliptic equation similar to (|3.4p also arises in the 
study of the self-similar solutions of the scmilinear parabolic equation 

« t - Am = in M™ x M+. 



See the work of Peletier, Terman, and Weissler [2lJ and Haraux and Weissler [141 ] . 

The theory of principal eigenvalues of fully nonlinear operators goes back to the 
work of Lions [l|| , who used stochastic methods to study the Bellman equation in a 
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bounded domain. Recently, several authors, including Birindelli and Demengel [H, 
[f| , Ishii and Yoshimura fig , and Quaas and Sirakov p2| have studied the principal 
eigenvalues of more general fully nonlinear operators in bounded domains (see also 
[l|). Our methods in this section are similar to those in these works, although 
extra complications arise from the unboundedness of R™. The special form of (|3.4p . 
particularly the gradient term, provides us with enough compactness to overcome 
these obstacles. 

R, r > 0, and assume that u G C(M. n \B r ) satifies 

VI A {D 2 u) - • Du < au inW l \B r , 
u < on dB r , 

and u(y) < Ce~ a ^ P for some constants C, a > and p > 1 . Then there exists a 
constant R = R(a, A) > such that r > R implies that u < 0. 

Proof. Set R := 2(a + A + 1). A routine calculation verifies that the function 
ip(y) := exp(— |y|) satisfies 

(3.6) P- (D 2 ^(y)) - X -y ■ Dcp(y) > (-A + \y\/2) y{y) > (a + l)<p in {|y| > R}. 

Suppose on the contrary that there exists \yo\ > r such that u(yo) > 0. Let 

T] := inf {p > : pip > u in « n \B r } . 

Then ij > 0, rjip > u in R n \B r , and owing to the faster decay rate of u, there exists 
| j/i | > r such that rjip(yi) — u(yi). In particular, the map 

y i — ► u(y) — nip{y) has a local maximum at y = y\. 

Recalling that u is a viscosity solution of (|3.5p . we see that 

P~ (r/L»V(yi)) - V^Vi ■ D( p{yi) < au{yi) = rjaifiyr) < r]{a + l)<p(yi), 



a contradiction to (J3T6J) . □ 

The next lemma is an analogue of [J Theorem 3.3], adapted to our setting. 
The result goes back to observations due by Berestycki, Nirenberg and Varadhan 
Q, and earlier versions have been developed by Quaas and Sirakov e 22], and Ishii 
and Yoshimura [l6| to study the principal half-eigenvalues of fully nonlinear elliptic 
operators in bounded domains. It is an important comparison tool which is essential 
to our analysis in this section. 

Lemma 3.2. Assume that u, v, f S C(R™) ; / > 0, and ael satisfy 

F(D 2 u) -^y-Du-au<f< F(D 2 v) - -y ■ Dv - av in R". 

Suppose also that v > and u(y) < Cexp(— a|y| p ) for some p > 1 and every 
y € R™. Then either u < v in R™, or u = tv for some constant t > 1. 

Proof. Let R — R{a 1 A) be as in Lemma 13. II For s > 1, define w s := u — sv. Then 
p- A (D 2 w s ) - -y ■ Dw s -aw s <f- sf < in R™. 
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For any s > 1 so large that w s < on Br , Lemma 13.11 implies that w s < on 
M. n \Bn. Hopf's lemma then implies that w s < in R". Define 

t := inf {s > 1 : w s < in R n } . 

Then t > 1 and wt < 0. If t = 1, then u < u in R™, and we have nothing left to 
show. Suppose that t > 1. We claim that w t = 0. If not, then by Hopf's lemma, 
wt < 0. Select 5 > so small that t — S > 1 and wt-s < on Br. Then wt_a < 
on R n , as we argued above. This contradicts the definition of t, completing the 
proof. □ 

Corollary 3.3. Suppose that u, v E C(M. n ) and a£l satisfy 

F(D 2 u) -^y-Du-au<0< F(D 2 v) - -y ■ Dv ~ av in R™. 

Suppose also that v > and u(y) < C exp(— a|j/| p ) for some p > 1 and every 
y G R". TTien either u < in R", or u = tv for some constant t > 0. 

Proof. Suppose that w(yo) > for some g/o € R". Using the homogeneity of F and 
by multiplying u by a large positive constant, if necessary, we may assume that 
u(yo) > v(yo). According to Lemma [IOJ there exists t > such that v = tu. □ 

In our analysis, an important role with be played by the functions {4> a } a>0 , 
which we define by 

4> a (y) :=exp(-a\y\ 2 ). 
Lemma 3.4. For a = 1/4A, the function <p a satisfies 

nX \ j.a / -n- ,n2j.a\ 1 ,. n^o/ /( n_1 ) A , 1 



(3-7) ^j^n-A^-^^^T- ^ 
Likewise, for a — 1 /4A, t/ie function <j) a satisfies 

(n-l)A l\ ,„ , . , 1 ,„ /rtA 



(3.8) + - j r < v« ~- 2 y°r< K 2A 

Proof. For any a > 0, the Hessian of a is given by 

D 2 r(y)=r(y) (4a 2 y(g>y-2al) 1 

and has eigenvalues (4a 2 |y| 2 — 2a)cf> a (y) with multiplicity 1 and —2acf> a (y) with 
multiplicity n — 1 . Hence 

rr> 2 <6 a r )i = J a ^W An ~ 4aA M 2 )> fOT a11 \v\ < ( 2a rK 

A ' Al (y>> \a<p a (y)(2\(n- 1) + 2A - 4aA|y| 2 ), for all |y| > (2a)"i 
Therefore, for each a > 0, 

(3.9) P A - A ( J D 2 a )-i 2/ .^ a = a (2An-4 a A|y| 2 + |y| 2 )^ for | v | < (2a)-*, 

and 
(3.10) 

P A - A (C 2 Q )-iy^ a = a (2A(n - 1) + 2A - 4aA|y| 2 + |y| 2 ) Q for | V | > (2a)"i 

In particular, for every a < (4A) , we have 

(3.11) V^ h {D 2 <t> a ) - -y ■ D<j> a > {2a\n)(j) a in R". 
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For a = (4A) we carefully check that 



(3.12) v^ K {p*<ir)--vD<if< , )A 

Similar calculations demonstrate that 

(3.13) V^ K (D 2 (j) a ) — • Z30 a = a (2An — 4aA|y| 2 + |y| 2 ) (j) a for < (2a)"*, 

and 
(3.14) 

p+ A (£> 2 ^)- iy.£>^ = a (2A(n - 1) + 2A - 4aA|j/| 2 + |y| 2 ) «T for \y\ > (2a)~*. 
Thus for any a > (4A) _1 , we have 

(3.15) ^a^V) - ^ • Dcf) a < (2aAn)(j) a in R" . 
Finally, for a = (4A) -1 , we calculate 

(3.16) v+ A (D 2 r)~lyDr> (^r^+l)^ inR "- D 

For the rest of this section, we let b denote the special constant 

From (|3.9|) and (|3.10[) , we see that 4> b satisfies 

(3.17) V-(D^ b ) \y ■ Dcf> b > + ^\V\ 2 ) / in 
Define a norm ||-||„ on C(R n ) by 

|M(y)|exp(6|y| 2 ) 

and let X denote the Banach space 



sup 

1 + M 2 



I = {!i£C(l") : Hull, <(»}. 
Notice that convergence in X implies uniform convergence in R n . Define the set 

C:= {ueX :0< u{y) < C exp(-%| 2 ) for some C > 0} . 
Notice that C is a convex subset of X. 

Proposition 3.5. For each v G X such that v > 0, there exists a unique solution 
u *E X of the equation 

(3.18) F(D 2 u) -\j ■ Du = v inW n . 

Moreover, u G C. 

Proof. We will first demonstrate existence. For each R > 0, let u R be the unique 
solution of the Dirichlet problem 

F (D 2 u R ) - • Du B . = v in B R , 
ur • = on dBn. 
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Set K := \\v\\^ max{4A/raA, 16A}. According to (13.17|) . the function tp :— K(j) b is a 
supersolution of 

F(I>V) -\yDi>> ||«||, (1 + |y| 2 ) <p b > v in M". 

By the maximum principle, < < ip = K e^cp(—b\y\ 2 ) fo r every i? > 0. Using 
local C Q estimates for uniformly elliptic equations (c.f. [261]), we deduce that for 
each fixed Rq > 0, 

Sup ||«Jj||cf«(B Ro ) < 00. 
R>R + 1 

Extend ur to be zero outside K"-\Br, and extract a subsequence Rj — ► oo such 
that 

— * w locally uniformly in R™ 

for some function u £ C(M"). Evidently < u < ip, and thus u £ C. From the 
stability properties of viscosity solutions under uniform convergence, it follows that 
it is a solution of equation (|3. 18|) . Uniqueness follows from Corollary 13.31 Indeed, 
if u%, «2 £ C are solutions of (|3.18|) . then the function w := U\ — U2 satisfies 

P~(D 2 w)- -y-Dw< inM™. 

Comparing w with <f> 2b and using Corollary [331 we deduce that w < in R™. □ 
We denote by X + the set 

X + = {u G X : u > 0} . 

Let „4 : AT + — > A" + be the solution operator of (|3.18[) . That is, A(v) :— u, where u 
is the unique solution of (|3.18p . Then A (X + ) C C. It will be convenient to use the 
notation 

F a [u] := F(D 2 u) -^y-Du- em. 

Define the constant 

(3.19) a+(T) := sup {a : there exists ip G A"+\{0} such that F a [<p] > in M™} . 

We call a + (-F) the positive anomalous exponent of i* 1 . When there is no ambiguity, 
we will drop the dependence on F and write a + = a + (F). From Corollary 13 . 31 and 
Hopf 's Lemma, it is clear that the anomalous exponent satisfies 

a + < inf {a : there exists ip £ X + \{0} such that F a [ip] < in W 1 } . 

From Lemma [32] we see that 
(3.20) 

nA _i_ \ (n-l)A + A n (n — 1)A + A j_,^4_ , nA 
2A * a+ faO * 2A * 2 ^ 2A * " +( ^ ) * 2A ' 
Moreover, if A ^ A, then all of the inequalities in (|3.20[) are strict. If F and G are 
two uniformly elliptic, positively homogeneous operators such that F < G, then 

a+(F) < a+{G). 

Thus from (|3.20|) we deduce that for every operator F satisfying (|2.2[) and (|2.3[) . 

nA I . „, nA 
(3.21) _< a+( ir)<_. 
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We will show that a + is an eigenvalue of the operator A, using the Leray- 
Schauder alternative. For the convenience of the reader, we first state this result. 
A proof can be found in [l3j |. 

Definition 3.6. If Y and Z are Banach spaces, we say a (possibly nonlinear) map 
A : Y — > Z is compact if, for each bounded subset B C Y, the closure of the set 
{A(x) : x € B} is compact in Z. 

Theorem 3.7 (Leray-Schauder Alternative). Suppose Y is a Banach space, and 
C C Y is a convex subset of Y such that S C. Assume that A : C — > C is a 
(possibly nonlinear) junction which is compact and continuous. Then at least one 
of the following holds: 

(i) the set {ieC : x = fiA(x) for some < /i < 1} is unbounded in Y , 

or 

(ii) i/iere exists x f= C /or which x — A(x). 

In order to apply Theorem 13. 71 in our setting, we must verify that the nonlinear 
operator „4 is continuous and compact. 

Proposition 3.8. The operator A is continuous and compact with respect to || • || ^ . 

Proof. Let {vk}k>i Q A such that UufcU, < 1. Set Uk ■= A{vk). Let e > be 
given, and fix a large constant R > to be selected below. As in the proof of 
Proposition I3.5[ from (|3.9p and (13.10|) we see that the function <fi := M<j> b satisfies 

V-{D 2 cf) -ly.Dcf>>(l + \y\ 2 ) 4> b > v k in K" 

for M := max {4A/nA, 16A} and every k > 1. It follows that u k < i/j for every 
k > 1. Using local C Q estimates, we have 

SUP |K || c°(B R ) < °°- 

fe>l 

Therefore, we may select a subsequence, which we also denote by k, such that 

lim sup ||tt fc - ui\\ l ^i Br ) = 0. 

Now take R = (2M/e) 1/2 . Then for any y > R and k, I > 1, 

Ml/) - ex P (&|y| 2 ) < 2|V»(y)|exp(6|^| 2 ) = 2M < g 



l + |y| 2 ~ 1 + i? 2 1 + i? 2 

It follows that 

lim sup ||ttfc - uiW^ < e. 

K^oo k t l>K 

A diagonalizing procedure now produces a subsequence of {uk} which is Cauchy in 
X. Therefore, A is compact. 

To see that A is continuous, suppose in addition that the sequence Vk converges 
strongly in A to a function ud. In particular, Vk — > v uniformly in R". We can 
find u € X and a subsequence such that Uk j — > u in A, and hence uniformly 
in R n . By the stability properties of viscosity solutions with respect to uniform 
convergence, it follows that u = A(v). By uniqueness, the full sequence Uk converges 
to u. □ 
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Proposition 3.9. There exists a unique ip + G X such that f + (0) = 1 and 

(3.22) F(D 2 Lp + ) - ■ D(p+ = a + Lp + m E™. 

Moreover, ip + G C n C^" (M™ ) , and <y9 + S Cj^"(R n ) i/F is concave or convex. 
Proof. Select u> £ C such that ||w||^ = 1. We claim that for any e > 0, 

(3.23) if it G C and a > satisfy u = aA(u + then a < a + . 
Indeed, for such a > and u G C, we have 

F Q [u] > in R". 

Since w ^ 0, if a > then it ^ 0. We see from definition (|3.19[) that in this case 
a < a + . Obviously if a = 0, then u = 0. Our claim (|3.23|) is confirmed. 
We may now apply Theorem 13.71 to deduce that for each e > 0, the set 
D e := {u G C : there exists < a < a + + e such that u = aA(u + sw)} 
is unbounded in X. Select u e G D e such that ||iie||* > 1- Let a e > such that 

u £ = a e A{u £ + ew) 

Evidently, a e > 0. Normalize by setting v £ := u e / Hwell*; and notice that by the 
homogeneity of A, the function v e satisfies 

v £ = a £ A(v £ + ew/ \\u £ \\J . 

By the compactness of A, we may select ip + G X, a number < a* < a + , and a 
subsequence Sj — > 0, such that 

v £j — > </? + in X and a ej — » a*. 

Since A is continuous, it follows that <p + = a*A(if + ). Thus cp + G C. Clearly 
||v + ||* — 1j an d thus a* > 0. By Hopf's Lemma, yj + > in R". 

We will now argue that a* — a + , and that (^ + is unique up to multiplication by 
a positive constant. Suppose that a > a* and ip G X + \{0} are such that 

F(D 2 -0) - • Dip > aii inR n . 

Then 

(3.24) F Q . [<^+] = < F a ty] < F a , ty] . 

According to Corollary 13.31 we have i\> = tip + for some t > 1. This implies that 
a* = a. Recalling (|3.19|) . we see that a* > a + . Recalling that by construction 
a* < u + , we deduce that a* = a + . Moreover, we have shown that tp + is unique 
up to multiplication by a positive constant. 

The last statement in the proposition follows from the standard regularity theory 
for uniformly elliptic equations (c.f. @, @|). □ 

Proof of Theorem Define 

(3.25) <f>+(x,t) := t~ a+ tp+ 

Assuming that $ + and ip + are smooth, and using (|3.22|) . we easily verify that $ + 
is a solution of (| 1 . 3|) . If $ + and ip + are not smooth, our calculation can be made 
rigorous in the viscosity sense by the use of smooth test functions. The uniqueness of 
a + and <I> + is established by performing this computation in reverse and appealing 
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to Proposition 13.91 All of the corresponding assertions regarding a and $ now 
follow from Remark I2JJ □ 



We conclude this section with an estimate of our self-similar solution <i> + from 
above and below, and an example. 

Lemma 3.10. For each < a < (4A) _1 , there exists a constant C > such that 

(3.26) ^+(y)<Cexp(-aM 2 ). 

Likewise, for each a > (4A) _1 , there exists a constant C > such that 

(3.27) exp(-a|y| 2 ) <C<p+(y). 

Proof. By construction, since <p + £ C we have that the estimate (|3.26j) holds for 
ai = (8A) _1 . We will therefore only show (|3.27[) . as a similar argument obtains 
(pT26f for all ai < (4A)" 1 . For a > (4A)~\ by (|3T4)) we have that 

V + {D 2 4> a ) - -y ■ D<j> a < a + 4> a in R n \B r , 
provided that we take r > so large that 

a(2A(n- 1) + 2A) - a+ 



r 2 > (2a)- 1 and r 2 > 



a(4aA- 1) 



Also take r > R, where R = R (a + , A) is the constant in Lemma [37X1 Let C be so 
large that cf> a < C(p + on B r . Then the function w := <fi a — Ctp + satisfies 

V^(D 2 w)-^y Dw <a + w in R n \B r , 

and w < on dB r . According to Lemma T3.ll the function w < in R™. That is, 

(f> a < Cip + in M" . □ 

Corollary 3.11. For each < a\ < (4A) _1 < (4A) _1 < a%, there exists a constant 
C > 1 such that 

(3.28) C-H- a+ exp(-a 2 \x\ 2 /t) < $+(aj, t) < Cr a+ exp (~ ai \x\ 2 /t) 

for all (x,t) E M™ x R+. 

Example 3.12. Consider the case that F is convex. Then F is a supremum of 
a collection of linear operators L k with constant coefficients, and each of which 
satisfy (|2.2[) and (|2.3[) . Since a + (L fc ) = a~(L k ) = n/2 for every fc, we deduce that 

77 

«-(F)<-<a+(^). 

We claim that these inequalities are strict unless F is linear. Suppose that a + (F) = 
n/2. Let <p and ifk be the functions obtained in Proposition 13.91 for F and L k , 
respectively. Notice that 

K/2 M=0 = L k n/2 [<p k ] < F n/2 [ip k ] in R". 

According to Corollary 13. 3[ tp = ifk for every k. That is, the fundamental solutions 
of the constant-coefficient linear parabolic operators L k are equal. This implies 
that L k — L for every k (see Friedman [IH). Hence F = L. 



ASYMPTOTICS FOR FULLY NONLINEAR PARABOLIC EQUATIONS 



13 



Remark 3.13. In the case that F(M) depends only on the eigenvalues of M, we 
immediately deduce that ip + and hence <&+(•,£) are radial functions. This follows 
from the invariance of the equation under an orthogonal change of variables, and 
Corollary |3.3l In particular, the self-similar solutions corresponding to the operators 
V + and V~ are radial. 



4. Asymptotic convergence to self-similar solutions 



In this section, we present the proof of Theorem 11.21 Owing to Remark l2.ll we 
need only prove the first statement. For ease of notation, we write a = a + (F) and 
$ = $ + . Fix a solution u = u(x, t) of the equation 

(4.1) ut + F(D 2 u) = in R" x R + , 
subject to the initial condition 

(4.2) u(x,0)=g(x). 

We require the initial data g to be continuous, not identically zero, and to satisfy 
the condition 

< g(x) < C e- B ^ 2 
for some constants B, Cq > 0. For a > 0, we denote 

u° := %u(x, t) := a a u (a 1/2 x, at) . 

For each a > 0, the function u a is a solution of (|4.1[) . 

We intend to show that as the parameter a — > oo, the rescaled solutions u a 
converge locally uniformly in R n x R + to a positive multiple of Q(x, t). Recall that 
$ is invariant under T a : 

(4.3) ®(x, t) = T CT $(x, t) = cr Q $ (v 1/2 x, at) for all (x, t) e R" x R + , a> 0. 



The proof of Theorem 11.21 will consist of a scries of lemmas. As a preliminary 
step, we show that u is bounded between positive multiples of possibly shifted 
in time. 

Lemma 4.1. For each t > there exists a constant C > 0, depending only on 
Cq and t, such that 

(4.4) u(x,t) < C$(x,t + T) for all (x, t) € R n x R+. 
Proof. If t > Tgv, then according to Corollarv l3.11l 

$(x,r) > c{r)e- B ^\ 

provided we choose c(r) > small enough. Thus C$(x, r) > g(x) for C := Cq/c{t) 
and x £ R™. The maximum principle implies that CQ(x,t + t) > u(x,t) for all 
(x,i) G R n x R + . □ 

Lemma 4.2. For each to,T > 0, there exists C > 0, depending only on B, Co, to, 
and t, such that 

(4.5) <£>(x,t) <Cu(x,t + r) for all x£W l ,t>t . 
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Proof. By the strong maximum principle, u(x, t) > on E™. Let C > be so large 
that 

t ) < Cu{x, t + t) for all |x| < 1, 

and 

$(x,t) < Cu(x,t + r) for all |jc| = 1, < t < to- 
Applying the maximum principle, we have <&(a;,io) < Cu{x,to + t) for all x G K™, 



and (|4.5|1 follows from another application of the maximum principle. □ 
For r = 1/(2XB), we use (I4.3[) to rewrite the inequality (I4.4[) in terms of vF as 

(4.6) u CT (:M) < C1>(x,t + r/a) for all (xi)ei"xR +1 oO. 

Recalling (|3 . 28[) . we see that for some constant C(t) > depending only on a lower 
bound for t > 0, in addition to B and Co, we have the estimate 

(4.7) u"(x,t) < C(t)exp(-|a:| 2 /8A(t + r/cr)) for all (x,t) G R n x K + , cr > 1. 

According to (|4.7[) and local Holder estimates for solutions of uniformly parabolic 
equations (see Wang 25|, Theorem 4.19]), we obtain 

SU P 1 1 u<7 1 1 C-<(Q) < 00 

(T>1 

for some < 7 < 1 and any compact parabolic domain Q C M™ x K + . Therefore, 
for every sequence Ok — > 00, we may select a function t/ € C(R" x E+) and a 
subsequence, also denoted by tJfe , such that u ah — » C/ locally uniformly in K" x M + . 
By the stability of viscosity solutions with respect to uniform convergence, each 
such rescaled limit U is a solution of equation (|4.ip . 

Let S denote the set of such sequential limits {U} of the family {u a } a >i. We 
will prove Theorem 11.21 by showing that S is a singleton set consisting only of a 
positive multiple of $. 

Lemma 4.3. There exists a positive constants C > such that for all U G S, 

(4.8) C" 1 ^ < U < C$ itil"xl + . 

Proof. From gU), we see that for each U G S and (x, t) G R" x R + , 
U{x,t) < limsupu <T (x,t) < C$(x,t) . 

a — >oo 

For the other direction, fix to,r > 0. According to (|4.5|) , for each t > the 
inequality 

(4.9) o- a <$> [a 1/2 x, at) < Ca- a u (a 1/2 x, at + r) 

holds for all (x, t) G W 1 x R + and all sufficiently large a > 1. Rewrite (|4.9[) as 

$(a;,t) < Cu' 7 + r/cr) . 

Take the lim-inf of the right side as a — > 00 to see that $(2:, t) < CU(x,t) for any 
U eS. □ 

Notice that Lemma |4~31 and local Holder estimates imply that 

(4.10) sup II^Hc^q) < 00, 

ues y ' 

for every compact subset Q C R n x M + . 
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Define the constant 

(4.11) C* := inf {C > : there exists U € S such that U < C<S>} . 

In light of Lemma l4~3l < C* < oo. We will eventually show that S = {C*$}. 
There are two basic steps in the proof. First, we will show that U < C*<I> for every 
U E S. Second, we show that if U ^ C*$ for some U E S, then we can find another 
function V E S and a small number S > such that V < (C* — <5)$, in contradiction 
to the definition (|4. 1 1|) of C* . Most of the subtlety in the proofs of these statements 
arise from difficulties in managing the "tails" of <&. These obstructions are removed 
by the construction of a special subsolution, which we use as a comparison function. 

Lemma 4.4. For each a > 4?, there exists r, r\ > and a subsolution w — w(y, s) 
of the differential inequality 

(4.12) w s +V^ A (D 2 w)-^yDw<0 mR n xK+, 
satisfying the initial conditions 

w(y,0) < e~ a \ y \ for all \y\ < r, and w(y,0) < — 7/e~' y for all \y\ > r, 
and such that for each R > there exists a time S > such that 
w(y, s) > for every \y\ < R, s > S. 

Proof. Select a constant a > (4A) _1 and set <p(y) := </>°(y) = e~ a \ y \ . Recall from 
(f3T5| that 

V + (D 2 ip) ~^y- Dtp< (2akn) <p in R". 

Let (3 := l+2aAn and r\ := 2(/3+A+l), and define a function ip{y) = min {e _ri , e - ^}. 
Recalling (|3 . 6[) and that the minimum of supersolutions is a supersolution in the 
viscosity sense, we see that 

V-(D 2 i>)-^y-Di)>0 mB ri+1 , 

and 

V~ (D 2 ijj) -^y-DiP>({3 + l)ip in R n \B ri . 
Now define (p(y,s) := e~ ,3s p(y) and ij(y, s) := —er^ +1 ' s tp{ij). Then 
(p s + V + {D 2 (p) - • Dtp < -<p inR"xK + , 

and 

S + V + (D 2 i>) -^y-D4><{(3+ l)e- {l3+1)s e- ri x in R" x R + , 

where x = 1 on B ri and % = on W l \B ri . Set 

1 



(3+1 

Then Cp{y, s) > 6(f3 _+ l) e - ( -> 3+1 ^- ri for all y E B ri and s > 0. Therefore, the 
function w := <p + Stp satisfies (|4. 12[) . 

We now investigate the set of (y, s) for which w > 0. For every \y\ > n, 



w(y, s) = e-P* (e~ a ^ 2 - 5e~ s e-^\ 
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From this expression, we observe that w(y, s) > whenever s > a\y\ 2 — \y\ + \ogS 
and \y\ > n. Finally, select r > large enough that r > r\ and as 2 > s + log § 
whenever s > r. This choice of r ensures that 

w(y,0) < --e - ' 1 ' 1 for all |y| > r. 

Taking rj := 5/2, the proof is complete. □ 

Corollary 4.5. There exist r,rj > swc/i i/iai /or any R > and any subsolution 
u of 

u t +V~{D 2 u) <0 ffll"x(l,oo) 
satisfying initial conditions 

u {x, 1) < —1 /or every |x| < r, and u(x, 1) < yye - ' 31 ' /or every \x\ > r, 
i/iere exists T > 1 swc/i i/iai 

t) < for all t > T, \x\ < R^ft. 
Proof. Let r, n, and w be as in Lemma T4. 41 for a = ^j. Define 

u(a;, t) := -tu f -^=,log*j , (x, ()6l"x (l,oo). 

Then v is a supersolution of the equation 

v t +V~(D 2 v) > in R" x (1, oo), 

such that 

v(x, 1) > ne~' x ' for every \x\ > r, and v(x, 1) > — e - "' 1 for every \x\ < r. 

In particular, v > u at time t = 1 and thus v > u in 1" x (1, oo) by the maximum 
principle. According to the conclusion of Lemma T4.4[ for each R > there exists 
S > such that 

w ( x ,t) = -w (t~ 1/2 x,logt) < for all r 1/2 \x\ < R, log t > S. 
Hence the conclusion is obtained for T = exp(S). □ 
Lemma 4.6. For any U G S, 

(4.13) U(x, t) < C*$(x, t) for all (x, t) E R n x E+. 

Proof. Let r, n be as in Corollary 33] and fix a small number e > 0. Recalling (|4.7[) . 
we may choose constants C\ , a > such that 

it CT (:r, 1) < Ci exp(-a|a;| 2 ) for all a; £ R", cr > 1. 

Set m := minui< r &(x, 1) and select r\ > r such that 

— e- a|a:|2 < 7ye- |K| for all Id > n. 
me 

According to the definition (|4.11|) of C*, we may select ci > 1 such that 

u ai (x, 1) < (C* +e)$(x,l) for all |x| < n. 
Define w := — (C* + 2e)<i>. Then w is a subsolution of the parabolic equation 
w t +V^{D 2 w) < inR"x(l,oo), 
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and at time t = 1 the function w satisfies 

w(x, 1) < — 1 ) < —me for all |x| < r, 
w(x, 1) < -e$(x, 1) < for all |x| < n., 

and 

w(x, 1) < u ai (x, 1) < (me)r)e~\ x \ for all |x| > r%. 

According to Corollary |4.51 for each R > there exists a time T = T(R) > 1 such 
that 

w(x, t) < for all |x| < itt 1/2 and t > T. 

This reads 

of u (cr^x, o-!^ < (C* + 2e)$(x, t) for all |x| < Rt 1 ' 2 and t > T. 
Thus for any a > o\ , 

u a {x,t) = a a u (a 1/2 x,at) 

= a a u (al /2 {a/a 1 ) 1 / 2 x,a 1 (a/a 1 )tj 

< {a/a 1 ) a (C* + 2s) $ (V/cxi) 1 /^, (a/ai)t) 
= (C* + 2e) $ 0, t) 

provided that (cr/cri) 1 / 2 ^ < R (at/ai) 1 ^ 2 and at/a± > T(R). In particular, for 
each (x,t) € R™ x R + , there exists 02 > <J\ large enough that 

it CT (x, t) < (C* + 2e) $ [x, t) for all a>a 2 - 

It follows that for any U € S, 

£/(x,<) < (C* + 2e)$(x,t) for all (i,i)el"xl + . 

The conclusion is obtained by sending e — * 0. □ 



To complete the proof of Theorem 11.21 we will need elementary properties of S 
contained in the following two lemmas. 

Lemma 4.7. IfUeS, then T a U 6 S for any a > 0. 

Proof. Select a sequence aj — * oo such that u aj — > U locally uniformly in R n x R + . 
It is easy to check that for aj := aaj, the sequence u aj — > T a U locally uniformly in 
R™ x R + . □ 

Lemma 4.8. The set S is closed in the topology of local uniform convergence. 

Proof. Let Uj £ S such that Uj — > U locally uniformly in R n x R + . Fix a compact 
subset Q of R™ x R + , and select aj — > oo such that 

sup|u CTj -Uj\< 2-i. 
Q 

It is clear that u aj converges to U as j ' — ► oo, uniformly on Q. Now a diagonalization 
argument produces a sequence dj — > oo for which the functions u aj converge to U 
locally uniformly in R™ x R + , as j — > oo. □ 

Lemma 4.9. Suppose that U E S and C > are such that U < C<I?. TTien either 
U = C$ or tfiere exists (5 > and e S such that V < (C — <$)<&. 
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Proof. Suppose that U £ S and C > are such that U < C$, but U ^ C$ in 
K" x M.+ . By the strong maximum principle, U(x,t) < C$(x,i) for every (x,t) £ 
K™ x M. + . Let r, rj > be as in Lemma 14.41 and choose e > so small that 

U(x, 1) < (C — e)$(x, 1) for every \x\ < r. 

Let m := min| x |< r 1). Select 6 > small enough that 

^— - > — and <5 1/2 $(x, 1) < e" |a;| for all b| > r. 
Denote by w the function 



which is a subsolution of the equation 

w t +V~(D 2 w) < 



Moreover, 
and 



w 



(x, 1) < (5 1 / 2 $(x, 1) < e ! x ! for every |x| > r, 



is 

w(x, 1) < T-pr— t) < — 77 for every lacl < r. 



According to Corollary I4.5[ for any R > there exists > 1 such that 

U(x, t) <{C- 5) $(x, t) provided that |x| < iiVt, and t > T. 

It follows that for each (x, t) £ W 1 x M+ there exists cr' > 1 large enough that 

%U{x,t) < (C-5)<f>(x,t) for all cr > cr'. 

According to Lemma B~T1 T CT (J7) e S. Recalling (pLTTOj) . we may select V G C(E™ x 
M+) such that up to a subsequence, 7^£7 — > locally uniformly in 1™ x K + . It is 
clear that 

V<(C-S)9. 

According to Lemma I4T51 V £ S. □ 

Proof of Theorem \1.SX According to Lemmas 14.61 and 14.91 and the definition (|4.11|) 
of the constant C* , the function C*$ is the only element of S. The proof of 
Theorem II .21 is complete. □ 



Remark 4.10. Let us repeat a remark made in |17j. If we express the constant 
C* obtained in Theorem 11.21 as a function of the initial data, C* — C* [g] , we 
see immediately that a nonnegative solution u of (|1.3p has the property that t i— > 
C*[u(-, t)] is constant. We thereby deduce a conservation law for our fully nonlinear 
equation, generalizing the conservation of mass in the case of a linear operator. 
It would be interesting to discover more information about C* [g] in the general 
nonlinear case. What is this conserved quantity? 
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